Large Fourier transforms never exactly realized by braiding conformal blocks 
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Fourier transform is an essential ingredient in Shor's factoring algorithm. In the standard quantum cir- 
cuit model with the gate set {U(2), CNOT}, the discrete Fourier transforms Fn = (lu 13 )n xn ,i, j = 
0, 1, ■ ■ ■ ,N — 1, u> = e~ , can be realized exactly by quantum circuits of size 0(n 2 ),n = logiV, and so 
can the discrete sine/cosine transforms. In topological quantum computing, the simplest universal topological 
quantum computer is based on the Fibonacci (2+l)-topological quantum field theory (TQFT), where the stan- 
dard quantum circuits are replaced by unitary transformations realized by braiding conformal blocks. We report 
here that the large Fourier transforms Fn and the discrete sine/cosine transforms can never be realized exactly 
by braiding conformal blocks for a fixed TQFT. It follows that approximation is unavoidable in implementation 
of the Fourier transforms by braiding conformal blocks. 



Introduction. The simplest topological model for quan- 
tum computing which can approximate any quantum circuit 
efficiently by braiding conformal blocks is based on the Fi- 
bonacci topological quantum field theory (TQFT) [3]. The 
corresponding conformal field theories (CFTs) for the Fi- 
bonacci TQFT include the level=l WZW G 2 CFT. TQFTs 
are low energy effective theories for topological phases of 
matter such as fractional quantum Hall (FQH) liquids, where 
quasi-particles can be anyons, even non-abelions theoretically. 
We will use the term anyon loosely here to include also non- 
abelions. On theoretical and numerical grounds it is believed 
that the Fibonnaci TQFT is an essential part of an effective 
theory for the FQH liquids at filling fraction v = 12/5 1 10] 
jj]. Moore and Read proposed that the ground state wave- 
functions for anyons localized at fixed positions are given by 
the conformal blocks of the corresponding conformal field 
theory |5]. Thus quantum gates in topological quantum com- 
puters are the braiding matrices of the conformal blocks, 
which are also the braiding statistics of anyons. 

A decade ago, Shor discovered the polynomial-time quan- 
tum algorithm for factoring integers. A key component of 
Shor's algorithm is the application of the discrete Fourier 
transforms Fn- It is known that the Fibonacci topological 
quantum computer can simulate Shor's algorithm efficiently, 
but the simulation requires approximations of the Fourier 
transforms 1 3]. In this paper we present a "no-go" theorem by 
showing that approximation is unavoidable. Closely related to 
the Fourier transforms are the discrete sine/cosine transforms 
which are also useful for signal processing. Our discussion 
for Fourier transforms applies equally to those transforms. 

As the prospect of a topological quantum computer has 
attracted increased attention, examination of the program- 
ming and compiling issues attendant to this design has begun 
la]. Even an accurate (10 -5 ) NOT gate requires several hun- 
dred elementary braids according to the known approximation 
scheme |9]. Audiences seeing such compilations always ask, 
"yes, but isn't there a better way? Can't the arithmetic proper- 
ties of Fibonacci anyons be matched to the number theory of 



factoring?" While efficient factoring is still a theoretical pos- 
sibility, we show no arithmetic wizardry will create the all- 
important Fourier, sine or cosine transforms inside TQFTs. 

A TQFT has a finite label set L = {a, b, c, • • • }, which 
physically represents the anyon types in the theory. Then a 
TQFT is a consistent rule to assign each 2-dimensional ori- 
ented compact space £ a vector space V(S), and each cobor- 
dism (M,Si,E 2 ) a linear map ^(M,S 1 ,S 2 ) : V(£i) — > 
V(Yj2). In particular, a projective representation of the map- 
ping class group M(£) on V(E). When E has boundaries, the 
boundaries will be labelled by anyons. 

A TQFT is unitary if each vector space V(E) has a posi- 
tive definite Hermitian inner product < •, • >s satisfying the 
following conditions: 

1) : The Hermitian inner product is multiplicative with re- 
spect to disjoint union of surfaces, and the inner product on 
F(0) for the empty surface is 1. 

2) : The Hermitian inner product is natural with respect to 
the mapping class group action. 

3) : For any cobordism (M, Ei,E2) and any x € V(Ei) 
and y S V(E 2 ), we have 

< Z(M, Ei, E 2 )(a;), y >e 2 = < x, Z(M, E 2) S x )(y) > Sl . 

Theses conditions imply that the projective representations 
of the mapping class groups are unitary. Furthermore, accord- 
ing to II ill for any TQFT and any surface E (if <9E ^ 0, 
then <9£ should be labelled) a spanning set for V(S) is ob- 
tained by the functor V applied to 3 -manifolds M contain- 
ing a labelled trivalent graph with dM — E. Thus for any 
x,y E 7(E) with x = Z(M), we have < Z(M),y > s =< 
Z(M,%,dM){l),y > E =< X,Z{M,dM,$)(y) > . It fol- 
lows from this identity that any Hermitian structure obey- 
ing l)-3) above is determined by the operators Z(M, dM, 0), 
hence unique. We can use the gluing axiom to reduce the com- 
putation of the Hermitian inner products for all surfaces to the 
computation for annuli and pairs of pants. It follows that if all 
the quantum dimensions of an Hermitian TQFT are positive, 
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FIG. 1: F-moves: the free ends are labelled by anyons 



and the Hermitian products on all pairs of pants are positive 
definite, then the TQFT is unitary. 

F-matrices Given a unitary TQFT and a 4-punctured 
sphere S^ bcd , where the 4 punctures are labelled by anyons 
of types a,b,c, d. The 4-punctured sphere can be divided 
into two pairs of pants(=3-punctured spheres) in two differ- 
ent ways. In FIG.l, the 4-punctured sphere is the boundary of 
a thickened neighborhood of the graph in either side, and the 
two graphs encode the two different pants-decompositions of 
the 4-punctured sphere. The F-move is just the change of the 
two pants-decompositions. 

By the axioms of a TQFT, each pants decomposition 
of Sl bcd determines an orthonormal basis of V(S% b c d ). 
Therefore the F-move gives rise to a change of orthonormal 
bases of the same Hilbert space V(S^ b c d ), hence induces a 
unitary matrix F a>b ^ d , which is called the F-matrix. 

From the definition the F-matrices are unitary, but it is not 
obvious that the entries of the F-matrices are always algebraic. 
One of our goals is to show that the entire unitary structure, 
including the F-matrices, is compatible with algebraic choices 
for all unitary TQFTs. The difficulty lies in the choices of the 
F-matrices as they are basis dependent. The obvious "solu- 
tion": solve all complex equations in a TQFT such as the pen- 
tagon and hexagon equations for their real and imagine parts 
independently plus the unitarity constraints for the F-matrices, 
is not sufficient for the compatibility as the condition of be- 
ing purely real or purely imaginary is not algebraic. Our ap- 
proach instead is to satisfy the algebraic conditions first for 
the F-matrices with certain normalization and then to deduce 
unitarity from the normalization. 

Fibonacci TQFT First we recall the data for the Fibonacci 
TQFT, our chief example. There is only one non-trivial anyon 
type t in the theory. We will also use r to denote the golden 
ratio r = , and no confusions should arise. 

There are two unitary TQFTs with anyon types {1, r} and 
the fusion rule: r ® r = 1 © r. One is the mirror (or parity 
reversed) theory of the other. We list the data for one theory 
and refer to the resulting theory as the Fibonacci TQFT. The 
data for the other theory is obtained by complex conjugate all 
the data below. 

Anyon types: {1,t} 

Fusion rule: 1 (g> r = r ® 1 = t,t <B> t — 1 © r 
Quantum dimensions: {l,r} 
Twists: #i = 1, 6 T = e~ 
Braidings: K[ T = e~ , i?^ r = e~ 
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FIG. 2: The basis is in one-one correspondence to admissible la- 
bellings of the internal edges with 1 or r subject to the fusion rules 
at each trivalent vertex. In all figures, label represents type 1, and 
label 1 represents r. 




1 



1 



1 



FIG. 3: The braiding is obtained by stacking the braid on top of a 
conformal block basis, and the braining matrix is computed by using 
the graphical calculus. 
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S-matrices: S\ = ^ 2+T I |..v= < 

Topological degeneracy: 

Let Yj g n be the genus=g oriented surface with n boundaries 
labelled by r, then dimV{Y. g , n ) = r " +( ~ 1) " 7 ^ 

Topological inner product: 

The Hilbert space V(T, gtn ) is spanned by labelled uni- 
trivalent graphs {G} in a bounding handlebody H g . n (for sim- 
plicity we ignore the framing subtlety.) Given two vectors 
in v, w G V(T,g t n) represented by two graphs G Vl G w , then 
the inner product of v, w is the topological invariant of the 
3 -manifold M with a trivalent graph G inside obtained from 
doubling the handlebodies and uni-trivalent graphs G V ,G W : 
glue the orientation reversed handlebody containing G v with 
the handlebody containing G w by the identity map on their 
boundaries. 

Conformal block basis (FIG. 2): 

r -l r -l/2N 

The braiding of two anyons in a conformal block basis state 
is represented by the following graph (FIG. 3): 

To find the matrix elements, we form the inner products of 
this braided basis with all basis. The topological inner product 
in the conformal block basis is given by flipping over the first 
argument and stacking on top of the second argument. Hence 
the matrix element is an invariant of a trivalent graph with cer- 
tain braidings. Now we observe that the invariant of any such 
graph is a complex number in the number field Q(y/r, £20). 
where £20 = e 27 ™/ 20 and Q(y / r, £20), an example of a num- 
ber field, consists of all complex numbers which are rational 



F-matrices: F = 
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polynomials in \fr, £20 with integer coefficients. Therefore, 
we have: 

Observation: All matrix entries of the braiding matri- 
ces with the above choices of data for the Fibonacci TQFT 
lie inside the number field Q(y/r,£2o) whose Galois group 
is the non-abelian dihedral group D4. Furthermore, only 
1,2,4,5, 10,20-th roots of unity exist in Q(-/r, £20). 

We will now see that there are only finitely many roots of 
unity in Q{\/t,£2o)- But to realize all the discrete Fourier 
transforms Fff, we need infinitely many root of unity, there- 
fore discrete Fourier transforms Fn for large N cannot be 
realized exactly by braiding conformal blocks. The roots of 
unity in £20) determine which Fourier transform can 

be potentially realized by braiding conformal blocks in the Fi- 
bonacci TQFT. The notation below is clarified in the follow- 
ing section. Notice that [Q(y/r,£ 2 o) : Q] = [Q(Vt,&o) = 
Q (£20) ][Q (£20) : Q]. For a primitive m-th root of unity £ m , 
[Q(£m) : Q] = <P{m), where <fi(m) is the Euler function 
whose value is the number of integers from 1 to m — 1 that 
is relatively prime to m. Since r £ Q(^2o)> so [Q(\Zt, £,20) ■ 
Q(6o)] = 2, and [Q(6o) : Q] = 0(20) = 8. Hence we have 
[QCv^i 60) : Q] = 16. It follows that there are only finitely 
many roots of unity in Q(y / r, £20) since there are only finitely 
m such that <fi(m) < 16. 

Another consequence of this observation is that the Fi- 
bonacci TQFT cannot be realized in an abelian extension of Q 
because the Galois group of Q(\/t, £20) is non-abelian. It suf- 
fices to show that the Galois group of Q(y/r, i) C Q(Vr, £20) 
is non-abelian, which is the same as the Galois group for the 
minimal polynomial of y/r f(x) = x 4 — x 2 — 1. To deter- 
mine the Galois group of f(x), we use the following fact: let 
g(x) — x 4 + ax 2 + b 6 Q(x) be irreducible with Galois group 
G. If neither b nor b(a 2 — 46) is a square in Q, then the Ga- 
lois group of g(x) is the non-abelian dihedral group D4. For 
a proof, see Proposition 4.1 1 of |6] on Page 273 and Ex. 9 on 
page 277. Now it is obvious that the Galois group of f(x) is 

Finally let us determine all the possible roots of unity 
in Q(\Zt,^2o)- If a primitive m-th root of unity £ m is in 
Q(\/t, £20), then 4>(m) is a factor of 16 because [Q(\/t, £20) : 
Q] = [Q(V?,Cao) : 0(£m)][Q(£m) : Q] = 16. If mis rela- 
tively prime to 20 and > 7, then £2o m would be in Q(\/t, £20)- 
But (7j(20m ) = 80(777) > 16, a contradiction. It follows that 
if £ m € Q(\/t, £20), then 777 is of the form 2 k ■ 3 • 5 for pos- 
sibly k = 1,2, 3. But first 3 cannot be a factor m because 
otherwise, £qo £ Q(y / T, £20)- Since Q(^6o) would be a sub- 
field of Q(a/t, £20) which are both degree 16 extension of Q, 
we will have Q(^6o) == QiV^ , £,20) ■ But this is impossible 
since the Galois group of Q(^6o) is abelian, while the Ga- 
lois group of Q(\/t, £20) is non-abelian. Exactly the same 
argument will rule out k — 3 with £40 replacing £qq. So 
the only possible primitive roots of unity in Q(\/t, £20) are 
£ m ,m — 1, 2, 4, 5, 10, 20, and their powers. 

Using the relation = e 27ri c ' 8 , we deduce that the cen- 
tral charges of the corresponding CFTs are c = 14/5 mod 8, 
which is realized by the level=l G2 CFT Because the cen- 



tral charges c ^ 0, we have to either work with projective 
representations rather than linear representations of the map- 
ping class groups or work with some central extension of the 
mapping class groups for extended surfaces. For the torus 
case, the projective representation can always be lifted to a lin- 
ear representation as follows: direct computation shows that 
(srj) 3 = £±-s 2 , so if we set t = t(^-)~^, then (st) 3 = s 2 . It 
has been shown that a 3rd-root of unity of is sufficient to 
lift all projective representations of the mapping class groups 
to linear representations of the extended mapping class groups 
full . Hence there are at least three different normalizations 
for a given TQFT, which lead to successively larger number 
fields: 

1) . Arbitrary choice for the F-matrices and projective rep- 
resentations for the mapping class groups 

2) . Unitary normalization for the F-matrices and projective 
representations for the mapping class groups 

3) . Unitary normalization for the F-matrices and linear rep- 
resentations for the extended mapping class groups. 

For the Fibonacci TQFT, with normalization 1), the Fi- 
bonacci TQFT can be described in Q(^2o); with normaliza- 
tion 2), Q(v^7 £20); with normalization 3) Q(\/t, £eo)- Note 
this field contains £ m for all m|60 by an argument similar to 
the one above. 

Unitary TQFTs Let Q be the field of the rational numbers 
(a field here is not in the sense field theory in physics, but as 
in number theory. A field is a generalization of the number 
systems Q, R, C.) A number field is a finitely dimensional 
vector space K over Q which is a field: a vector space with 
a compatible multiplication. The field K is called an exten- 
sion field of Q, and the dimension of K as a vector space over 
Q is called the degree of the extension, denoted by [K : Q] . 
Given a complex number x, Q(x) is the field of all complex 
numbers of the farmp(x)/q(x), where p(x), q(x) are polyno- 
mials in x with coefficients in Q and q(x) 7^ 0. For example, 
Q(\/t) is a degree=4 extension of Q. Fields can be extended 
repeatedly as follows: let K be an extension of Q and y a 
complex number, then K (y) is the field of all complex num- 
bers of the form p(y)/q(y), where p(y), q(y) are polynomials 
in y with coefficients in K and q(y) ^ 0. Given two com- 
plex numbers x, y, the number field Q(x, y) is the extension 
of first Q to Q(x) = K or Q(y) = K, then K to K(y) or 
K(x), which both are Q(x,y). The degree of the extension 
is [Q(x,y) : Q] = [Q(x,y) : Q(x)][Q(x) : Q] = [Q(x,y) : 

Q(y)][Q(y):Q]. 

Main Results: 1. Given a unitary TQFT, there is a nor- 
malization so that all the entries of the F-matrices are in a 
number field K, and the F-matrices associated to the F-moves 
are unitary. 

2. Each Hilbert space V(E) has an orthonormal basis so 
that every representation matrix of the mapping class group 
has entries in the number field K. Warning the Galois group 
Gal(K/Q) is not necessarily abelian. 

3. Large Fourier transforms, the discrete sine/cosine trans- 
forms cannot be realized exactly in any fixed TQFT by braid- 
ing conformal blocks. 
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Parts 3 follows from Part 2 as follows. We recall that the 
number of roots of unity in a number field K is always finite. 
To see this, the degree of the extension of Q(£ m ), £ m = e"^ 
is 4>(m), where <f){m) is the Euler function. Therefore if [K : 
Q] = n and <f>(m) > n, then £ m cannot be in K because 
otherwise the extension degree will be > n. Similarly for 
Q(sm(27r/m)),Q(cos(27r/m)). 

Part 2 follows from Part 1 . Given a unitary TQFT, there is 
a unique way to construct compatible topological inner prod- 
ucts for V(S)'s r fiUl . Chapter IV, Section 10], and we need 
an explicit orthonormal basis for each V(H) to compute the 
braiding matrices. To do this one sets up a graphical calculus 
so that each matrix entry is an invariant of a certain trivalent 
graph, which depends on our choices of the F-matrices and 
9 symbols. The theorem is reduced to the careful choices of 
F-matrices and 9 symbols which are compatible with the topo- 
logical inner product. The invariants of such graphs are poly- 
nomial of certain roots of unity and 6j symbols F^ n . There 
are three kinds of contributions of roots of unity: the braiding 
eigenvalues, the twists, and the higher Frobenius-Schur indi- 
cators resulting from bending anyon trajectories. They are in 
some fixed extension of Q whose degree is determined by the 
fusion rules through Vafa's theorem. The 6j symbols are con- 
strained by the pentagon identities. To have a consistent set 
of 6j symbols F^ n with graphical calculus, it is sufficient to 
solve the following set of polynomial equations (for easiness 
of notation we drop the dependence on trivalent vertices): 

i pijk _ I _dM_s:. 

o pijm _ jpklrn" _ pmij I d„, d n 
^■ r kln — r ijn" ~ r nk'l*\J djdi 

J - 2-~in kp*n mns* r lkr* q"kr* mis* 

Any solution of this set of equations will be a consistent 
choice of 6j symbols for the unitary TQFT. Now we cite a 
theorem in algebraic geometry: the solution to the polyno- 
mial equations above is an algebraic variety over Q(y/di),i = 
1 , 2, • • ■ , R, where R is the number of anyon types. Since this 
variety has at least one point which gives rise to the TQFT, 
then there will be also an algebraic point by Theorem 7 on 
Page 32 1 7]. It follows that every graph invariant will be in- 
side a fixed finite extension of Q(^fdi) and hence in a number 
field over Q. 

The resulted graphical calculus from the solution of the 
the pentagon equations with the above normalization has very 
nice properties. The conformal block basis is an orthogonal 
basis. The 9 symbols 6(a, b, c) = y/d a dbd c , where d a , db, d c 
are the quantum dimensions of the anyons a,b,c. One con- 
sequence of the 9 symbol values is that the conformal block 
basis elements have the same length, independent of the inter- 
nal labellings. So the ^-matrices are change of basis for two 
orthonormal bases up to overall scalars, hence are unitary. 

Approximation by Fibonacci Quantum Computer Since 
the exact realization of the Fourier transforms is impossible in 
the Fibonacci TQFT, we would like to approximate them us- 
ing braiding matrices. Given a prescribed accuracy, it will be 
interesting to find the explicit approximations. We will only 
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1 0/1 10/11 1 0/1 1 0/1 1 

FIG. 4: The edges labelled by or 1 correspond to standard qubits, 
and other basis span the non-computational subspace which should 
be evolved by the identity operator ideally in the computational pro- 
cess. Again represents 1, and 1 represents r 

outline an approximation here. 

To simulate a standard n-qubit quantum circuit Ul ■ 
(C 2 )®" -> (C 2 )®", we embed (C 2 )®" into the conformal 
blocks on In + 2 Fibonacci anyons at fixed positions. Since 
dim(V2n+2) = F > 2n+2 > 2" except for n = 1, we need to 
choose an efficiently computable subspace of the conformal 
blocks. One way to do this is to choose the following sub- 
space (C 2 )®" of V~2n+2 with the conformal block basis (FIG. 
4). 

Then we look for a braid b so that the following diagram 
commutes up to the prescribed error, where p(b) is the braid- 
ing matrix of conformal blocks. 

(C 2 )®" — > V 2n +2 
Ul I I P {b) 

(C 2 )®" — » V 2n+ 2 

The standard quantum circuits for the exact realization of 
the Fourier transforms are given on Page 219 of [9]. Given 
a precision e > 0, then one finds a braid that approximate 
Fjy by using the approximations of the single qubit gates and 
CNOT in |1. 

Conclusion TQFTs are effective theories for topological 
phases of matter such as the fractional quantum Hall liquids. 
Specifically, the braiding matrices of conformal blocks are 
unitary transformations of the degenerate ground states when 
anyons are fixed at certain positions. Because polynomial 
time approximation schemes exist |9], the reported obstruc- 
tion to exact realization of the Fourier transforms will not im- 
pose a fundamental physical constraint on topological quan- 
tum computing. However as a practical matter there is an im- 
portant distinction between billions as opposed to millions of 
braid generators to factor a large number. 

The Jones braid representation(s) that we get from Fi- 
bonacci anyons can be described as a regularized Fourier 
transform "FTB" of the braid group(s) B n . The braid gen- 
erators correspond to "position" coordinates and the path ba- 
sis of conformal blocks is a regularized momentum basis for 
the group algebra of the braid group C[B n ]. The chosen reg- 
ularization consists of passing to an appropriate semi-simple 
quotient, the Temperley-Lieb algebra TL^ 1 = C[B n }/ ~, q = 
e 2iri/5 j^yg snown t na j one cannot find the FT of large 
cyclic groups inside these FTB. The most direct application 
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of FTB is to the estimation of Jones polynomial evaluations 
J3l and 01 ■ The possibility of harnessing FTB for number 
theoretic application such as factoring should be explored. 
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